LONG-TIME ASYMPTOTICS FOR THE KORTEWEG -DE VRIES 
EQUATION VIA NONLINEAR STEEPEST DESCENT 



KATRIN GRUNERT AND GERALD TESCHL 



Abstract. We apply the method of nonlinear steepest descent to compute the 
long-time asymptotics of the Korteweg— de Vries equation for decaying initial 
data in the soliton and similarity region. This paper can be viewed as an 
expository introduction to this method. 



1. Introduction 

One of the most famous examples of completely integrable wave equations is the 
Korteweg-de Vries (KdV) equation 

(1.1) q t (x,t) =6q(x,t)q x (x,t) - q xxx (x,t), (i,i)elxl, 

where, as usual, the subscripts denote the differentiation with respect to the corre- 
sponding variables. 

Following the seminal work of Gardner, Green, Kruskal, and Miura [T7J, one 
can use the inverse scattering transform to establish existence and uniqueness of 
(real-valued) classical solutions for the corresponding initial value problem with 
rapidly decaying initial conditions. We refer to, for instance, the monographs by 
Marchenko [57] or Eckhaus and Van Harten [T^ . Our concern here are the long-time 
asymptotics of such solutions. The classical result is that an arbitrary short-range 
solution of the above type will eventually split into a number of solitons travelling 
to the right plus a decaying radiation part travelling to the left, as illustrated in 
Figure [TJ The first numerical evidence for such a behaviour was found by Zabusky 
and Kruskal [39j . The first mathematical results were given by Ablowitz and Newell 
[I] , Manakov [55] , and Sabat [3D] . First rigorous results for the KdV equation were 
proved by Sabat [30] and Tanaka [34] (see also Eckhaus and Schuur [15], where 
more detailed error bounds are given). Precise asymptotics for the radiation part 
were first formally derived by Zakharov and Manakov [38 , by Ablowitz and Segur 
[5], [32], by Buslaev [5J (see also [S]), and later on rigorously justified and extended 
to all orders by Buslaev and Sukhanov 7 . A detailed rigorous proof (not requiring 
any a priori information on the asymptotic form of the solution) was given by Deift 
and Zhou [TUJ based on earlier work of Manakov [55] and Its [TO] (see also [50] , [5T] , 
[55]). For further information on the history of this problem we refer to the survey 
by Deift, Its, and Zhou [15]. 

To describe the asymptotics in more detail, we recall the well-known fact (see 
e.g. [9], [27]) that q(x,t) is uniquely determined by the (right) scattering data of 
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Figure 1. Numerically computed solution q(x, t) of the KdV 
equation at time t = 5, with initial condition q(x, 0) = sech(a; + 
3) — 5sech(a; — 1). 



the associated Schrodinger operator 



(1.2) 



H(t) 



dx 2 



+ q(x,t). 



The scattering data consist of the (right) reflection coefficient R(k, t), a finite num- 
ber of (t independent) eigenvalues —k 2 with < k\ < K2 < ■ ■ ■ < kn> and norming 
constants ~fj(t). We will write R(k) = R(k,0) and jj = jj(0) for the scattering 
data of the initial condition. Then the long-time asymptotics can be described by 
distinguishing the following main regions: 

(i). The soliton region, x/t > C for some C > 0, in which the solution is asymp- 
totically given by a sum of one-soliton solutions 



(1.3) 



q(x,t)~-2j2- 



^ 7 



~^ cosh (kjx — 4k jt — pj)' 
where the phase shifts are given by 



(1.4) 



Pj = \ log 



JV 



£ n 



J 1=3+1 



Kl 



In the case of a pure 7V-soliton solution (i.e., R(k, t) = 0) this was first established 
independently by Hirota [H] , Tanaka [33] , and Wadati and Toda [3S] . The general 
case was first established by Sabat [3U] and by Tanaka [33] (see also [TS] and [3T]). 

(ii). The self-similar region, \x/ (3t) x / 3 1 < C for some C > 0, in which the solution 
is connected with the Painleve II transcendent. This was first established by Segur 
and Ablowitz l32l. 



(iii). The collisionless shock region, x < and C 



< 



< C, for 



(3t)l/3(log(t))2/3 

some C > 1, which only occurs in the generic case (i.e., when R(Q) = — 1). Here 
the asymptotics can be given in terms of elliptic functions as was pointed out by 
Segur and Ablowitz [35] with further extensions in Deift, Venakides, and Zhou [T3j . 
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(vi). The similarity region, x/t < — C for some C > 0, where 

(1.5) q(x, t) ~ f^M^j ' sin (16ifc 3 - u(k ) log(192t*g) + S(k )), 
with 

^ )=-^log(l-|i?(fc )| 2 ), 

N 

*(*°) = 7 - arg(i?(fc )) + arg(r(ii/(fe ))) +4 Varctan(^-) 

3 = 1 

+ - [° log(|C-fc |)dlog(l-|i?(C)| 2 ). 

Here fc = a/ - t% denotes the stationary phase point, i2(fe) = R(k,t = 0) the 
reflection coefficient, and T the Gamma function. 

Again this was found by Zakharov and Manakov [38] and (without a precise 
expression for 6(ko) and assuming absence of solitons) by Ablowitz and Segur [2] 
with further extensions by Buslaev and Sukhanov 7 as discussed before. 

Our aim here is to use the nonlinear steepest descent method for oscillatory 
Riemann-Hilbert problems from Deift and Zhou [10] and apply it to rigorously 
establish the long-time asymptotics in the soliton and similarity regions (Theo- 
rem 231 respectively. 15.41 below). In fact, our main goal is to give a complete and 
expository introduction to this method. In addition to providing a streamlined and 
simplified approach, the following items will be different in comparison with [10j . 
First of all, in the mKdV case considered in [10] there were no solitons present. 
We will add them using the ideas from Deift, Kamvissis, Kriecherbauer, and Zhou 
[14] following Kriiger and Teschl [24] . However, in the presence of solitons there is 
a subtle nonuniqueness issue for the involved Riemann-Hilbert problems (see e.g. 
[H Chap. 38]). We will rectify this by imposing an additional symmetry condition 
and prove that this indeed restores uniqueness. Secondly, in the mKdV case the 
reflection coefficient R(k) has always modulus strictly less than one. In the KdV 
case this is generically not true and hence terms of the form R(k)/(1 — \R(k)\ 2 ) 
will become singular and cannot be approximated by analytic functions in the sup 
norm. We will show how to avoid these terms by using left and right (instead of just 
right) scattering data on different parts of the jump contour. Consequently it will 
be sufficient to approximate the left and right reflection coefficients. Details can be 
found in Section [6] Moreover, we obtain precise relations between the error terms 
and the decay of the initial conditions improving the error estimates obtained in 
Schuur [21] (which are stated in terms of smoothness and decay properties of R(k) 
and its derivatives). 

Overall we closely follow the recent review article [25] , where Kriiger and Teschl 
applied these methods to compute the long-time asymptotics for the Toda lattice. 

For a general result which applies in the case where R(k) has modulus strictly less 
than one and no solitons are present we refer to Varzugin [35] and for another recent 
generalization of the nonlinear steepest descent method to McLaughlin and Miller 
[28] . An alternate approach based on the asymptotic theory of pseudodifferential 
operators was given by Budylin and Buslaev [5] . 

Finally, note that if q(x,t) solves the KdV equation, then so does q(— x, — t). 
Therefore it suffices to investigate the case t — > oo. 
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2. The Inverse scattering transform and the Riemann-Hilbert 

PROBLEM 

In this section we want to derive the Riemann-Hilbert problem for the KdV 
equation from scattering theory. This is essentially classical (compare, e.g., [4]) 
except for two points. The eigenvalues will be added by appropriate pole conditions 
which are then turned into jumps following Deift, Kamvissis, Kriecherbauer, and 
Zhou [14j . We will impose an additional symmetry conditions to ensure uniqueness 
later on following Kriiger and Teschl [24 . 

For the necessary results from scattering theory respectively the inverse scat- 
tering transform for the KdV equation we refer to [27] (see also [4] and [9]). We 
consider real- valued classical solutions q(x,t) of the KdV equation (jl.ip . which 
decay rapidly, that is 



(2.1) max / (l + \x\)\q(x,t)\dx < oo, for all T > 0. 

Existence of such solutions can for example be established via the inverse scattering 
transform if one assumes (cf. [27l Sect. 4.2]) that the initial condition satisfies 

(2.2) f (l + \x\)(\q(x,0)\ + \q x (x,0)\ + \q xx (x,0)\ + \q xxx (x,0)\)dx<oo. 

JR 

Associated with q(x, t) is a self-adjoint Schrodinger operator 

d 2 

(2.3) H(t) = ——T+q(.,t), Q(H) = H 2 (R) C L 2 (K). 

dx z 

Here L 2 (R) denotes the Hilbert space of square integrable (complex-valued) func- 
tions over M and H k (R) the corresponding Sobolev spaces. 

By our assumption (|2.1| the spectrum of H consists of an absolutely continu- 
ous part [0,oo) plus a finite number of eigenvalues —k 2 e (— oo,0), 1 < j ' < N. 
In addition, there exist two Jost solutions ip±(k, x,t) which solve the differential 
equation 

(2.4) H{t)^±(k,x,t) = k 2 ^±(k,x,t), Im(fc) > 0, 
and asymptotically look like the free solutions 

(2.5) lim e Tikx *lj±(k,x,t) = 1. 

x— >±oo 

Both ip±(k,x,t) are analytic for Im(fc) > and continuous for Im(fc) > 0. 
The asymptotics of the two Jost solutions are 

(2.6) ^ ± ( k , x , t ) = e ±^(i + Q ± ( X]t )^L+0(^) 
as k — > oo with Im(A:) > 0, where 

(2.7) Q+(x,t) = - q{y,t)dy, Q_(x,t) = - q(y,t)dy. 

J x J —oo 

Furthermore, one has the scattering relations 

(2.8) T(k)ip T (k,x,t) = ij)±{k,x,t) + R±(k,t)ip±(k,x,t), k e K, 

where T{k) 1 R±(k,t) are the transmission respectively reflection coefficients. They 
have the following well-known properties: 
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Lemma 2.1. The transmission coefficient T (k) is meromorphic /or Im(fc) > with 
simple poles at 'iki, . . . ,ikn and is continuous up to the real line. The residues of 
T(fc) are given by 

(2.9) Res iKj T(k) = i^(t) 7+J (i) 2 = 
where 

(2.10) 7 + J W" 1 = ll^+K-,-,t)l| 2 

and ip + (\Kj , x,t) = fij(t)ip^(inj,x,t). 
Moreover, 



(2.11) T(k)R+(k,t)+T(k)R-(k,t) =0, \T(k)\ z + \R ± {k,t)\ z = 1. 

In particular one reflection coefficient, say R(k.t) = R + (k,t), and one set of 
norming constants, say "fj(t) = 7 +j(t), suffices. Moreover, the time dependence is 
given by: 

Lemma 2.2. The time evolutions of the quantities R(k,t) andjj(t) are given by 

(2.12) R(k,t) = R(k)e m3t , 

(2.13) lj (t)= lj e^\ 
where R(k) = i?(fc,0) and 7 j = 7 j(0). 

We will set up a Riemann-Hilbert problem as follows: 

(9 ^A\ m (k r t\ - / (r(fc)^-(fc,x,t)e ifc - V+(fc,^,t)e- ifca; ) , Im(fc) > 0, 
1 > 1 ' ' >~ \ (4>+(-k,x,t)e ikx T(-k)iJj_(-k,x,t)e- ikx ) , Im(fc) < 0. 

We are interested in the jump condition of m(k, x, t) on the real axis R (oriented 
from negative to positive). To formulate our jump condition we use the following 
convention: When representing functions on R, the lower subscript denotes the non- 
tangential limit from different sides. By m+(/c) we denote the limit from above and 
by m_ (k) the one from below. Using the notation above implicitly assumes that 
these limits exist in the sense that m(k) extends to a continuous function on the 
real axis. In general, for an oriented contour S, m + (k) (resp. m_(fc)) will denote 
the limit of m(/t) as k — > k from the positive (resp. negative) side of E. Here the 
positive (resp. negative) side is the one which lies to the left (resp. right) as one 
traverses the contour in the direction of the orientation. 

Theorem 2.3. Let S + (H(0)) = {R{k), k > 0; (Kj,7j), 1 < 3 < N} be the right 
scattering data of the operator H(0). Then m(k) = m(k,x,t) defined in (|2.14[) is a 
solution of the following vector Riemann-Hilbert problem. 

Find a function m{k) which is meromorphic away from the real axis with simple 
poles at Hnj and satisfies: 

(i) The jump condition 

, 01 . m ,. WM , M (l-\R{k)\ 2 -ip)e-**( fe )\ 

(2.15) m+(k) = m-{k)v{k), v(k) = I R ,Q e t*(k) \ )> 

for k G R, 
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(ii) the pole conditions 



Res iKj m{k) = Yon m{k) ^^(i^) q 

(2.16) 



fO _i 7 2g**(iKj) 

Res_ iK , m(k) = lim m(k) '3 

K—t—i-Kj \U U 

(iii) the symmetry condition 

1 



(2.17) m(-fc) = m{k) , ( () 

(iv) and i/ie normalization 

(2.18) lim m(i/c) = (1 1). 
_ffere £/ie phase is given by 

(2.19) $(fc) = 8ifc 3 + 2ifc-. 

Proof. The jump condition (|2.15|) is a simple calculation using the scattering rela- 
tions (|2.8p plus (I2.11|) . The pole conditions follow since T(k) is meromorphic for 
Im(fc) > with simple poles at ikj and residues given by (|2.9p . The symmetry 
condition holds by construction and the normalization (|2. 18|) is immediate from 
the following lemma below. □ 

Observe that the pole condition at \Kj is sufficient since the one at — 'iKj follows 
by symmetry 
Moreover, using 

(2.20) T(A0V-(fe, x, t)il>+(k, x,t) = l+ + (1.) 

as fc — > oo with Im(fc) > (observe that the right-hand side is just the diagonal 
Green's functions of Hit) divided by the free one) we obtain from (|2.6[) 

Lemma 2.4. The function m(k,x,t) defined in (|2.14p satisfies 

(2.21) m(k,x,t) = (l l)+Q( x ,t)^(-l 1)+0(^1 

Here Q(x,t) = Q+(x,t) is defined in (|2.7[) . 

For our further analysis it will be convenient to rewrite the pole condition as 
a jump condition and hence turn our meromorphic Riemann-Hilbert problem into 
a holomorphic Riemann-Hilbert problem following [T?]. Choose s so small that 
the discs |fc — 'iKj\ < s lie inside the upper half plane and do not intersect. Then 
redefine m(k) in a neighborhood of iKj respectively — iKj according to 



(2.22) m(k) 
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Note that for Im(fc) < we redefined m{k) such that it respects our symmetry 
(|2.17| . Then a straightforward calculation using Resi K m(k) — Hmk^i K (k — iK)m(k) 
shows: 

Lemma 2.5. Suppose m{k) is redefined as in (|2.22|) . Then m(k) is holomorphic 
away from the real axis and the small circles around ikj and —iKj- Furthermore 
it satisfies (|2. 15(1 . (|2.17|) . (|2.18p and the pole condition is replaced by the jump 
condition 



m + (fc)=m_(fc) ryV"''"^ , \k — in 



(2.23) 




(fc)=m_(fc) k+in, I \k + iKj\=e, 



where the small circle around 'iKj is oriented counterclockwise and the one around 
—iKj is oriented clockwise. 



Next we turn to uniqueness of the solution of this vector Ricmann Hilbcrt prob- 
lem. This will also explain the reason for our symmetry condition. We begin 
by observing that if there is a point fci € C, such that m(fei) = (0 0), then 
n (k) — fc _ 1 fci ro(fc) satisfies the same jump and pole conditions as m(k). However, 
it will clearly violate the symmetry condition! Hence, without the symmetry con- 
dition, the solution of our vector Riemann-Hilbert problem will not be unique in 
such a situation. Moreover, a look at the one-soliton solution verifies that this case 
indeed can happen. 

Lemma 2.6 (One-soliton solution). Suppose there is only one eigenvalue and 
that the reflection coefficient vanishes, that is, S+(H(t)) = {R(k,t) = 0, k € 
K; (n,j(t)), k > 0,7 > 0}. Then the unique solution of the Riemann-Hilbert prob- 
lem (|2.15j) - (|2.18|> is 9 iven by 
(2.24) mo(fc) = (f(k) f(-k)) 

I ^ + i K \-1..2 t$(i. 



^ ^ l + (2,c)-i 7 W> V + ^ (2K) 7 6 

Furthermore, the zero solution is the only solution of the corresponding vanishing 
problem where the normalization is replaced by lim^oo m(ifc) = (0 0). 
In particular, 

2 7 2 e **(i«) 



(2.25) Q(x,t) 



l + (2K)- 1 7 2 e**( iK )' 



Proof. By assumption the reflection coefficient vanishes and so the jump along the 
real axis disappears. Therefore and by the symmetry condition, we know that 
the solution is of the form mo(k) = (f{k) /(— fe)J where f(k) is meromorphic. 
Furthermore the function f{k) has only a simple pole at in, so that we can make 
the ansatz f(k) = C+D Z^f 1 . Then the constants C and D are uniquely determined 
by the pole conditions and the normalization. □ 

In fact, observe f(ki) = f{-k x ) = if and only if kt = and 2k = 7 2 e t<I, ( iK ). 
Furthermore, even in the general case m(ki) = (0 0) can only occur at k\ = as 
the following lemma shows. 
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Lemma 2.7. Ifm(ki) = (0 0) form defined as in (|2.14j) . then k\ — 0. Moreover, 
the zero of at least one component is simple in this case. 

Proof. By (|2.14[) the condition m(fci) = (0 0) implies that the Jost solutions 
tp-(k,x) and ij)+(k,x) are linearly dependent or that the transmission coefficient 
T{k\) = 0. This can only happen, at the band edge, k\ — or at an eigenvalue 
k\ = iKj. 

We begin with the case k\ — iKj. In this case the derivative of the Wronskian 
W(k) = ip+(k,x)ij/_(k,x) — ip' + (k, x)ijj- (k, x) does not vanish by the well-known 
formula j^W(k)\ k=kl = — 2ki Lip + (ki, x)tp_(ki, x)dx ^ 0. Moreover, the diagonal 
Green's function g(z 7 x) = W(k)~ 1 ip + (k, x)ip- (k. x) is Herglotz as a function of 
z = —k 2 and hence can have at most a simple zero at z = —k\. Since z — > —k 2 
is conformal away from z = the same is true as a function of k. Hence, if 
ij)+{\K,j, x) — ip—(iKj,x) — 0, both can have at most a simple zero at k = iKj. But 
T(k) has a simple pole at \Kj and hence T{k)ip_(k,x) cannot vanish at k = iKj, a 
contradiction. 

It remains to show that one zero is simple in the case fci = 0. In fact, one can show 
that 4rW{k)\k=ki ^ in this case as follows: First of all note that tp±(k) (where the 
dot denotes the derivative with respect to k) again solves Hip±(ki) = —k 2 ip±(ki) 
if fci = 0. Moreover, by W(ki) = we have VM-(fci) = cip-(ki) for some constant c 
(independent of x). Thus we can compute 

W(k 1 ) = W(ij; + (k 1 ),^(k 1 )) + W^ + (k 1 ),^(k 1 )) 

= c^W^+ih), ijj+ih)) + cW^-ih), ^_(fci)) 

by letting x — » +oo for the first and x — * — oo for the second Wronskian (in which 
case we can replace ip±(k) by e ), which gives 

W(h) = -iic + c- 1 ). 

Hence the Wronskian has a simple zero. But if both functions had more than simple 
zeros, so would the Wronskian, a contradiction. □ 

3. A UNIQUENESS RESULT FOR SYMMETRIC VECTOR RlEMANN— HlLBERT 

PROBLEMS 

In this section we want to investigate uniqueness for the holomorphic vector 
Riemann-Hilbert problem 

m+(fc) = m_(fc)w(fc), k G X, 

(3.1) m(— k) = m(k) 

lim m(in) = (l l) , 

K — KX> 

where we assume 

Hypothesis 3.1. Let E consist of a finite number of smooth oriented curves in 
C such that the distance between £ and {iy\y > yo} is positive for some yo > 0. 
Assume that the contour £ is invariant under k i— > — k and v(k) is symmetric 

(3.2) «("*)=(? J) w ( fc ) _1 (? J)' fceS - 
Moreover, suppose det(v(k)) = 1. 
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Now we are ready to show that the symmetry condition in fact guarantees unique- 
ness. 

Theorem 3.2. Assume Hypothesis \3.1l Suppose there exists a solution m(k) of 
the Riemann-Hilbert problem (|3.1[) for which m(k) = (0 0) can happen at most 
for k — in which case limsup fe _, TfTJFj * s bounded from any direction for j ' = 1 
or j — 2. 

Then the Riemann-Hilbert problem (|3.1|l with norming condition replaced by 
(3.3) lim m(in) = (a a) 

K — >00 

for given a£C, has a unique solution m a (k) — am(k). 

Proof. Let m a (k) be a solution of (|3.ip normalized according to (|3.3p . Then we 
can construct a matrix valued solution via M = (m, m a ) and there are two possible 
cases: Either det M(k) is nonzero for some k or it vanishes identically. 

We start with the first case. Since the determinant of our Riemann-Hilbert prob- 
lem has no jump and is bounded at infinity, it is constant. But taking determinants 
in 

M(-k) = M(k)(° 1 j). 

gives a contradiction. 

It remains to investigate the case where det(M) = 0. In this case we have 
W a (fc) = 5(k)m(k) with a scalar function <5. Moreover, S(k) must be holomorphic 
for k G C\S and continuous for k € £ except possibly at the points where m(ki) = 
(0 0) . Since it has no jump across S, 

5+(k)m+(k) = m Qj+ (fc) = m Qj _(fc)w(fc) = <5_(fc)m_(fc)v(fc) = <5_(/c)m+(fc), 

it is even holomorphic in C\{0} with at most a simple pole at k = 0. Hence it must 
be of the form 

5{k) = A + ?-. 

k 

Since 8 has to be symmetric, i5(fc) = 6(—k), we obtain B = 0. Now, by the 
normalization we obtain d(k) = A = a. This finishes the proof. □ 

Furthermore, the requirements cannot be relaxed to allow (e.g.) second order 
zeros in stead of simple zeros. In fact, if m(k) is a solution for which both com- 
ponents vanish of second order at, say, k = 0, then rh(k) = p-m(fc) is a nontrivial 
symmetric solution of the vanishing problem (i.e. for a = 0). 

By Lemma 12.71 we have 

Corollary 3.3. The solution m(k) = m(k,x,t) found in Theorem \2.3\ is the only 
solution of the vector Riemann-Hilbert problem (|2.15j) — 12. 18[> . 

Observe that there is nothing special about k — * oo where we normalize, any 
other point would do as well. However, observe that for the one-soliton solution 
(|2~24|) . f(k) vanishes at 

_ 1 - (2 K ) 2 7 V <I, ( iK ) 

" _ + (2 K )2 7 2gt<I>(i K ) 

and hence the Riemann-Hilbert problem normalized at this point has a nontrivial 
solution for a = and hence, by our uniqueness result, no solution for a = 1. This 
shows that uniqueness and existence are connected, a fact which is not surprising 
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since our Riemann-Hilbcrt problem is equivalent to a singular integral equation 
which is Frcdholm of index zero (see Appendix [Aj . 

4. Conjugation and Deformation 

This section demonstrates how to conjugate our Riemann-Hilbert problem and 
how to deform our jump contour, such that the jumps will be exponentially close 
to the identity away from the stationary phase points. Throughout this and the 
following section, we will assume that the R(k) has an analytic extension to a small 
neighborhood of the real axis. This is for example the case if we assume that our 
solution is exponentially decaying. In Section [6] we will show how to remove this 
assumption. 

For easy reference we note the following result: 

Lemma 4.1 (Conjugation). Assume that ECS. Let D be a matrix of the form 

where d : C\E — * C is a sectionally analytic function. Set 

(4.2) fh(k) = m(k)D{k), 
then the jump matrix transforms according to 

(4.3) v(k) = J D_(fc)- 1 u(fc)L>+(fc). 

If d satisfies d(—k) — d(k)^ 1 and lim^^oo d(in) = 1, then the transformation 
m(k) = m(k)D(k) respects our symmetry, that is, rh(k) satisfies (|2.17p if and 
only if m(k) does, and our normalization condition. 

In particular, we obtain 
respectively 

V 2 1 d + CL -dZ v 22 J 

In order to remove the poles there are two cases to distinguish. If Re(<5(iKj-)) < 0, 
then the corresponding jump is exponentially close to the identity as t — > oo and 
there is nothing to do. Otherwise we use conjugation to turn the jumps into one 
with exponentially decaying off-diagonal entries, again following Deift, Kamvissis, 
Kriecherbauer, and Zhou [T3]. It turns out that we will have to handle the poles 
at inj and — iKj in one step in order to preserve symmetry and in order to not add 
additional poles elsewhere. 

Lemma 4.2. Assume that the Riemann-Hilbert problem for m has jump conditions 
near \k and —in given by 

( 1 o\ 

m + (fc)=m-(fc)l i7 2 1 J, \k-iK\=£, 

m + (k) — m_(fc) I k +' 1K I , \k + ik\ = e. 
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Then this Riemann-Hilbert problem is equivalent to a Riemann-Hilbert problem for 
fh which has jump conditions near in and —in given by 



I (fc+i«) 

1 -TTJTT— 



(4.7) 



m+(k) = TO_(fc) ^ i7 3 (fe-ire) j ] 
( 1 °\ 

\ i 7 2(fe+i K ) L J 

d all remaining data conjugated (as in Lemma \4-1\ ) by 

k — \K Q 



Ifc — ikI = e, 



ik| = £, 



Proof. To turn 7 2 into 7 2 , introduce D by 



£>(fc) = < 



1 

ill 

N A:— iK 



k-\-iK 

k—'iK, 
k-\-in 





17-* 


1 



k-\-iK 
k — ik . 



k—'iK 
k-\-\K 






k— 



k-\-iK, 
k-i 



\k- 
\k- 

else. 



IK < £, 



IK < £, 



and note that D(k) is analytic away from the two circles. Now set m(k) = 
m(k)D(k)and note that D(k) is also symmetric. Therefore the jump conditions 
can be verified by straightforward calculations and Lemma 14.11 □ 

The jump along the real axis is of oscillatory type and our aim is to apply a 
contour deformation following [lOj such that all jumps will be moved into regions 
where the oscillatory terms will decay exponentially. Since the jump matrix v 
contains both exp(t$) and exp(— £$) we need to separate them in order to be able 
to move them to different regions of the complex plane. 

We recall that the phase of the associated Riemann-Hilbert problem is given by 



(4.8) 



$(Jfe) = 8ifc 3 + 2ik- 



and the stationary phase points, $'(fc) = 0, are denoted by ±fc , where 



(4.9) 



x 

12? 



For j > we have ko £ iR, and for | < we have fco € 
need the value kq defined via Re(<I>(iKo)) = 0, that is, 



For j > we will also 



(4.10) 



KQ 



> 0. 



We will set Ko = if j < for notational convenience. A simple analysis shows 
that for j > we have < k /i < k . 

As mentioned above we will need the following factorizations of the jump condi- 
tion (f2~T5|) : 



(4.11) 



v(k) = b-(k)- l b+(k), 
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where 



for |fc| > Re(fc ) and 

(4.13) «(*) = B_(fc)- 1 ( ri - | f* )|a ? W(fc), 
where 

/ i o\ A 

(4.14) B_(k) = g(fc)e^w 1 , = i-lW . 
for \k\ < Re(fco). 

To get rid of the diagonal part in the factorization corresponding to \k\ < Re(fco) 
and to conjugate the jumps near the eigenvalues we need the partial transmission 
coefficient T(k,ko). 

We define the partial transmission coefficient with respect to kg by 

k+iK i k Q e iR+ 

(4.15) T(k,k ) = { - N \ ( ko 




k 



for k <G C\S(fc ), where T,(k ) = [— Re(ko), Re(fco)] (oriented from left to right). 
Thus T(k, ko) is meromorphic for k € C\£(fco). Note that T(k, ko) can be computed 
in terms of the scattering data since |T(fc)| 2 = 1 — \R + (k, t)\ 2 . Moreover, we set 



(4.16) Ti(fco) 



Thus 



J2 2 «i> k o e «+, 

KjG(ko,oo) 
AT 

E + ^ L feo log(|T(C)| 2 R, ko E R+. 



(4.17) T(fc,fo) = l + Ti +o(JL), asfc^oo. 

Theorem 4.3. TTie partial transmission coefficient T(k, ko) is meromorphic in 
C\S(fc ), w/iere 

(4.18) E(fco) = [-Re(fc ),Re(fc )], 

with simple poles at iKj and simple zeros at —iKj for all j with Ko < Kj, and satisfies 
the jump condition 

(4.19) T + {k,k ) = T_(k,k ){l-\R(k)\ 2 ), forfceS(fco). 
Moreover, 

(i) T(-k, k ) = T(k, kp) - 1 , k g C\E(*o), 

(ii) T(—k, ko) = T(k, ko), k G C, m particular T(k, ko) is real for k e iR, and 
(hi) if ko e M+ i/ie behaviour near k = is given by T(k, k ) = T(k)(C + o(l)) 

witft C ^0 for Im(fc) > 0. 
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Proof. That \Kj are simple poles and — i/Cj are simple zeros is obvious from the 
Blaschke factors and that T(k,ko) has the given jump follows from Plemelj's for- 
mulas, (i), (ii), and (iii) are straightforward to check. □ 

Now we are ready to perform our conjugation step. Introduce 



1 



A>(*0, 



lKj\ < E, K < Kj, 



(4.20) D(k) = < 





k-\-'lK,j 



k-\-lKj 

1 



l-Do(fe), 
where 

'Tfcko)- 1 

T(k,k 

Observe that D{k) respects our symmetry, 



D a (k), \k + lKj\ < E, K < Kj, 

else, 



D (k) = 



D(-k) 



1 

1 



Now we conjugate our problem using D(k) and set 

(4.21) m(k) = m(k)D(k). 

Note that even though D(k) might be singular at k = (if fco > and i?(0) = —1), 
fh(k) is nonsingular since the possible singular behaviour of T(k, fco) -1 from Do(k) 
cancels with T(k) in m{k) by virtue of Theorem l4.3l (iii). 

Then using Lemma 14.11 and Lemma 14.21 the jump corresponding to kq < Kj (if 
any) is given by 

v(k) = P " 

(4.22) K - 

m= [ 

\ i 7 |e t ' , ' (i ' t j ) T(fe,feo) 2 

and corresponding to kq > Kj (if any) by 



; 2 e t*(.« J ) T(fcifeo) -2 



1 



1 

k+'iKj 



\k — iKj I 

|fc + iKj| = 



(4.23) 



S(ife) 



v{k) 



1 



i 7 2 e '*(f^T(fc,fco) 



\k — \Kj 



1 



In particular, all jumps corresponding to poles, except for possibly one if Kj = Ko, 
are exponentially close to the identity for t — * oo. In the latter case we will keep 
the pole condition for Kj — kq which now reads 



Res 



(4.24) 



Res_i Kj . m(fc) = lim m(k) 



-i7?e t *( i ^)T(ifi; : ,-,A;o)" 
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k e iK k Q e 





Figure 2. Sign of Re(<I>(fc)) for different values of k 



Furthermore, the jump along R is given by 
(4.25) v(k) 
where 



&_(A0 _1 &+(*)> fceK\E(fc ), 

s_(fc)- 1 s+(fc), fces(fc ), 



1 



(4.26) 6_ (fe) = ( q T ^Y°^ 



t 1 0^ 

M fc ) = fl(fc)c'*W 1 



and 



B-(k) = T_(fc,fc )- 



l-|H(fe)|- 



-7?(/c)e* <I, ( fe ) 1 



1 



■t$(fe) 




1 -^S)^(-fc)e- 4 * (fc) 



Here we have used 



R[-k) = R(k), k e 



T ± (-k,k )=T^(k,k )-\ fceS(fco) 



and the jump condition (|4.19p for the partial transmission coefficient T(k, ko) along 
S(fco) in the last step. This also shows that the matrix entries are bounded for 
k e M near k = since T±(-k, k ) = T±(k, k ). 

Since we have assumed that R{k) has an analytic continuation to a neighborhood 
of the real axis, we can now deform the jump along R to move the oscillatory terms 
into regions where they are decaying. According to Figure [5] there are two cases to 
distinguish: 

Case 1: fc € iM, fc ^ 0: 

We set £± = {k £ C| Im(fc) = ±e} for some small e such that £± lies in the 
region with ±Re(<f>(fc)) < and such that the circles around ±i«j lie outside the 
region in between £_ and £+ (see Figure [3]). Then we can split our jump by 
redefining fh(k) according to 



.27) 



m = 



mi^b+ik)- 1 , 

m(Jfe)&-(*)~S 
rh(k), 



< Im(fc) < e, 
—s < Im(fc) < 0, 
else. 
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Re($) > 



Re($) < 



Re($) > 



E_ 



/ Rc($) < \ 



Figure 3. Deformed contour for k e iK + 



Rc($)>0 



Ro(*)<0 



si \ 



Ro(*)<0 

si / si 



/— feo 



si / 

Re(*)>0_.'' 



Rc($)<0 



Si 

\, Rc(*)>0 



Figure 4. Deformed contour for fc S 



Thus the jump along the real axis disappears and the jump along E± is given by 

k e S+ 



(4.28) 



v 7 H-^)- 1 , k e E_ 



All other jumps are unchanged. Note that the resulting Ricmann-Hilbert problem 
still satisfies our symmetry condition (|2.17|) . since we have 



(4.29) 



By construction the jump along S± is exponentially close to the identity as t — ► oo. 
Case 2: fc G K, fc ^ 0: 

We set E± = E± U Ej_ according to Figure [4] chosen such that the circles around 
iiKj lie outside the region in between E_ and S + . Again note that Ej. respectively 
Ej. lie in the region with ± Re($(fc)) < 0. Then we can split our jump by redefining 
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i(k) according to 



(4.30) 



rh{k) 



m(k)b + (k)~ 1 , k between R and Si., 

m(fc)&_(fc) _1 , k between R and £]_, 

rh(k)B + (k)~ 1 , k between R and £+, 

m(fc)B_(fc)- 1 , /c between R and S?_, 

m(fc), else. 

One checks that the jump along R disappears and the jump along S-t is given 



by 

(4.31) 



v(k) 



'b+(k), ke^ + , 

b-ik)' 1 , fee Si, 

B+(k), k e S 2 + , 

.iL(fc)- 1 , fceiP. 



All other jumps are unchanged. Again the resulting Ricmann-Hilbert problem 
still satisfies our symmetry condition (|2.17j) and the jump along £±\{/co, — fco} is 
exponentially decreasing as t — > oo 



Theorem 4.4. Assume 
(4.32) 



(l + \x\) 1+l \q(x,0)\dx < 



for some integer I > 1 and abbreviate by Cj = i/ie velocity of the j 'th soliton 
determined by Re($(iKj)) = 0. Then the asymptotics in the soliton region, x/t > C 
for some C > 0, are as follows: 

Let e > sufficiently small such that the intervals [cj — e, Cj +e], 1 < J < iV, 
are disjoint and lie inside R + . 

If \j — Cj\ < e for some j, one has 



(4.33) 

respectively 
(4.34) 
where 
(4.35) 



q(v,t)dy 



q(x,t) 



27?(z.t) 



4 .^. Kl l + (2 Kl )-i 



i=j+l 



+ (2/ % )" 1 7?M) 



+ o(t- ; ), 



-4K/y?(x,t) 



{l + {2K j )-^Ax,t)f 



N 



0(t- 1 ), 



7 2(m)= 7 k 2 ^ +8 ^ n . 

4-+i Vrti+ ^ 



Iflj - Cj\ >£, for all j, 



has 



(4.36) 

respectively 
(4.37) 



J q(y,t)dy = ~4 J2 «i + 0(* ! ). Ko = \/J' 



(/«0 



g(a;,i) =0(t- ; ). 
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Proof. Since m(k) — fh(k) for k sufficiently far away from R equations (|2.2f I) . 
(|4.21| . and (|4.17|) imply the following asymptotics 

(4.38) m(fc) = (l l) + (-2T 1 (fc )+Q(x,f))^(-l l)+o(J-\. 

By construction, the jump along ~E± is exponentially decreasing as t — > oo. Hence 
we can apply Theorem IA. 61 as follows: 

^ If ~ c j\ > £ ( res P- l K o — K jl > £ ) f° r au J we can choose 7* = and 
u>* = w in Theorem \KM Since w is exponentially small as t — > oo, the solutions 
of the associated Riemann-Hilbert problems only differ by 0(t~ l ) for any I > 1. 
Comparing mo = (l l) with the above asymptotics shows Q+(x,t) = 2Xi(fco) + 
0(t ')• 

If || — Cj| < e (resp. |k§ — k 2 | < e) for some j, we choose 7* = 7^(2;, i) and 



to* = {y in Theorem I A. 6 ( where 



^(x, t) = 7 |e t *^)T(w i , i^=)- 2 = 7 | e - 2re ^+ 8 4* J] 

v3 ... 



w / x 2 

i >%- 



Ki + 

As before we conclude that w is exponentially small and so the associated solutions 
of the Riemann-Hilbert problems only differ by 0(t~ l ). From Lemma [2~6l we have 
the one-soliton solution mo(fc) = (f{k) /(— fc)) with 

l + (2/Cj) l ^(x,t) y k — iKj J 



As before, comparing with the above asymptotics shows 

27? (M) 
l+(2« i )-i 7 ?(x,t) 



q(x, t) = 2T X (fc ) + - ; , o ; ^ + o(*- 



To see the second part just use (|2.20p in place of (|2.21[) . This finishes the proof in 
the case where R(k) has an analytic extensions. We will remove this assumption in 
Section [5] thereby completing the proof. □ 

Since the one-soliton solution is exponentially decaying away from its minimum, 
we also obtain the form stated in the introduction: 

Corollary 4.5. Assume (|4.32|) . then the asymptotic in the soliton region, x/t > C 
for some C > 0, is given by 

N K 2 

(4.39) ,(*, t) = -2j2 CQsh2 , _ J 4k3 . } + 0(t~ l ), 

where 
(4.40) 








( Ki 









18 



K. GRUNERT AND G. TESCHL 




5. Reduction to a Riemann-Hilbert problem on a small cross 

In the previous section we have seen that for ko £ R we can reduce everything to 
a Riemann-Hilbert problem for fh(k) such that the jumps are exponentially close 
to the identity except in small neighborhoods of the stationary phase points ko and 
— ko. Hence we need to continue our investigation of this case in this section. 

Denote by £ c (±fc ) the parts of £ + U£_ inside a small neighborhood of ±fco. We 
will now show that solving the two problems on the small crosses £ c (fco) respectively 
£ c (— ko) will lead us to the solution of our original problem. 

In fact, the solution of both crosses can be reduced to the following model prob- 
lem: Introduce the cross £ = £1 U • • ■ U £4 (see Figure [5]) by 

Ei = {ue-™ /4 , u e [0, 00)} £ 2 = {uc i7r/4 , u € [0, 00)} 

(5.1) £3 = {we 3iTr / 4 , u G [0, 00)} £4 = {ue- 3 ™/ 4 , u e [0, 00)}. 

Orient £ such that the real part of k increases in the positive direction. Denote 
by D = {C, |CI < 1} the open unit disc. Throughout this section Q v will denote 
the function e lul ° e ^\ where the branch cut of the logarithm is chosen along the 
negative real axis (— 00, 0). 

Introduce the following jump matrices (vj for £ € £ 3 ) 

_ (l -i?i(C)C 2i,/ e~ t<E,(c) \ / 1 

Vl ~ {0 1 J ' V2 ~ \R 2 {OC Tl ^ mo 1 

( , _v _ (l -^(OC^e-^CON _ / 1 

(5-2) «3 " ^ 1 J, « 4 - (^ 4 ( C)C -2i V *(0 1 

and consider the RHP given by 

(5.3) M + (C)=M_(CK(C), CeS,, j = 1,2,3,4, 

M (C) -> I, C - 00. 

The solution is given in the following theorem of Deift and Zhou [TO] (for a proof 
of the version stated below see Kriiger and Teschl [25]). 

Theorem 5.1 ([TO])- Assume there is some po > such that Vj(() = I for |(j > po- 

Moreover, suppose that within |£| < po the following estimates hold: 
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(i) The phase satisfies $(0) = i$ G iK, $'(0) = 0, $"(0) = i and 

(5.4) ± Re (<i > ( C))>i|C| 2 , {+ l CeElUS3 ' 

(5.5) |a, (c) _$ ( o)-^|<qci 3 . 

(ii) There is some r£D and constants (a, L) € (0, 1] x (0, oo) such that Rj, 
j = 1, . . . , 4, satisfy Holder conditions of the form 

\Ri(Q-7\<L\{\", |ik(C)-r|<L|C| a , 

(5.6) |J? 3 (0 - Y^r-pl < £|Cl a , 1^(0 - Y^ypl < i|Cl a - 

TTzen £/ie solution of the RHP (|5.3p satisfies 



1 + a , 



/or | C| > po> where 

(5.8) /3 = v /^ e iW4-arg(r)+arg(r(i ! /))) e -ii# or i«/ ! y = _ J_ ^ /j _ 

Furthermore, if Rj(C) and ^(C) depend on some parameter, the error term is uni- 
form with respect to this parameter as long as r remains within a compact subset 
of D and the constants in the above estimates can be chosen independent of the 
parameters. 

Theorem 5.2 (Decoupling). Consider the Riemann-Hilbert problem 

m+(fe) = m-(k)v(k), k S S, 

(5.9) lim m(i«) = (1 1) , 

wii/i det(u) 7^ and fe£ < a < f3 < 2a, p(t) — > 00 &e given. 

Suppose that for every sufficiently small e > 6o£/i £/ie L 2 and £/ie L°° norms of 
v are 0{t~@) away from some e neighborhoods of some points kj £ S, 1 < j < n. 
Moreover, suppose that the solution of the matrix problem with jump v{k) restricted 
to the e neighborhood of kj has a solution which satisfies 

(5.10) M 3 (k)^l+-^^- + 0(p(t)-^, \k-k J \>e. 
Then the solution m(k) is given by 

(5.11) m(fc) = (l + ^E^+W)^), 
where the error term depends on the distance of k toT,. 

Proof. In this proof we will use the theory developed in Appendix lAl with mo(fc) = 
(l l) and the usual Cauchy kernel Q oa (s,k) = 1^^. Moreover, since symmetry 
is not important, we will consider C w on L 2 (Yj) rather than restricting it to the 
symmetric subspace L 2 (£). Here w± = ±(b± — I) correspond to some factorization 
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v = b_ 1 b+ of v (e.g., b- — I and b+ —v). Assume that m{k) exists, then the same 
arguments as in the appendix show that 

m(fc) = (l + — n(s)(w+(s) +w-(s))n ao (s,k), 
2tti Jx 

where p solves 

(l-C w )(p~(l l))=C ro (l 1). 

Introduce rh(k) by 

'mityMjik)- 1 , \k-kj\<2s, 
m(k), else. 



(5.12) m(k) 



(5.13) v(k) = 



The Riemann-Hilbert problem for m(k) has jumps given by 

'Mj(fc)- 1 , |k-fy|=2e, 
MjityvifyMjik)- 1 , k e E,e < < 2e, 

fc e s, |fc - %| < e, 

w(fc), else. 

By assumption the jumps are I + 0(p(t)~ a ) on the circles \k — kj\ — 2e and even 
I + 0(p{t)^ 13 ) on the rest (both in the L 2 and L°° norms). In particular, we infer 
that (I — C^) -1 exists for sufficiently large t and using the Neumann series to 
estimate (p - (l l)) = (I - Cy^C^ (l l) (cf. the proof of Theorem [OJ we 
obtain 



(5-14) l)ll^i3fc=W)--). 



C I|W||2 

clli&ll 

Thus we conclude 
(5.15) 



i(fc) = (l 1) + ^ / A(*)t»(«)' 



rf.s 



= (1 l) + ±J2f Ks)(M j (s)- 1 -I)^-+0(p(t)-^ 

Z7T1 j=1 J\s-kj\=2s S — K 

= (i i) +„<«)-» (i i)E^7+ w*)^). 

and hence the claim is proven. □ 

Now let us turn to the solution of the problem on S c (fco) = (£+ U £_) n 
{fc| |fc — fco| < e} for some small e > 0. Without loss we can also deform our 
contour slightly such that £ c (fco) consists of two straight lines. Next, note 

$(fco) = -16ifcg, $"(fc ) = 48ifc . 

As a first step we make a change of coordinates 

(5.16) C = y/48h(k-ka), k = k Q + —^= 



such that the phase reads $(fc) = <l>(fco) + h( 2 + 0(( 3 ). 
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Next we need the behavior of our jump matrix near ko, that is, the behavior of 
T(k, ko) near ko- 

Lemma 5.3. Let kg 6 R, then 

( k - kn \ w - 

(5.17) T(k,k )= (jjyj T(k,k ), 

where v = — — log(|T(fco)|) > and the branch cut of the logarithm is chosen along 
the negative real axis. Here 

N , . 

( , 18) fcM^-n^-p^, - vmwivc-*- 

is Holder continuous of any exponent less than 1 at the stationary phase point 
k = kg and satisfies \T(ko,ko)\ = 1. 

Proof. First of all observe that 

<" 9 > «* (si C MWMI 2 )^) - ' 

Holder continuity of any exponent less than 1 is well-known (cf. [29]). □ 
If k(() is defined as in (|5. 16[> and < a < 1, then there is an L > such that 

(5.20) T{k(C),k )-Cf(k ,k Q )e- i " lo ^ 2ko ^ S ^^ <£|Cf. 

where the branch cut of is chosen along the negative real axis. 
We also have 

(5.21) |#(C))- R(h)\ <L\C\ a 
and thus the assumptions of Theorem 15.11 are satisfied with 

(5.22) r = R{k Q )f(k , fc )-2 c 2i^log(2fcov / 48^) 

and v = — ^Mog(l — |i?(z )| 2 ) since |r| = |i?(z )|. Therefore we can conclude that 
the solution on Yf(ko) is given by 

(5.23) s v 7 



where /3 is given by 

/3 = ^ e i ^/ 4 ~ ars ^ +ars ( r ( i ^ ) )e~ t *( fco ^~ il ' 
{ "' 24) = V^c i ^/ 4 - arg ( fl ( feo » +arg ( r ^»¥(fc ,A : o) 2 (192fc3)-^e-*' I '( feo ^- i,y 

and 1/2 < a < 1. 
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We also need the solution M%{k) on S c (— fco). We make the following ansatz, 
which is inspired by the symmetry condition for the vector Riemann-Hilbert prob- 
lem, outside the two small crosses: 



(; ;)^c-*)(; j) 



(5-25) =1- - ) -4* ( ° n)+0(t- a ). 



Applying Theorem 15.21 yields the following result: 
Theorem 5.4. Assume 

(5.26) / (1 + |x|) 6 |g(a;,0)|dx < oo, 

Jr 

then the asymptotics in the similarity region, x/t < ~C for some C > 0, are given 
by 



/ q{y,t)dy = -A ]T ^ - - log(|T(C)| 2 R 



KjG(koi°o) 



(5.27) - y cos(16f*g - i/(*d) log(192i*g) + S(k )) + 0{r a ) 

respectively 

(5-28) q{x> t) = ^ ^(fco)fco gin(16tfc 3 _ ^ log(192ifco 3 ) + «5(fe )) + 0(r°) 
/or any 1/2 < a < 1. _ffere fco = ^ and 



(5.29) K*d)=--iog(|r(*o)|). 

7T 



5(fc ) -arg(i?(fc )) + arg(r(ii/(A:o))) +4Varctan(^ 



W- feo & \\T(ko)\ 2 J C-fco 
Proof. By Theorem 15.21 we have 



+ o(r a ), 

which leads to 

Q(x,t) = 2TM + -^-L(Re(/?)) + 0(t-«) 
upon comparison with (|4. 38[) . Using the fact that \fi/^/v\ — 1 proves the first claim. 
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To see the second part, as in the proof of Theorem 14. 4\ just use (|2.20j) in place 
of (|2.2ip . which shows 

q{x,t) = sJ^Im{(i) + 0{t- a ). 

This finishes the proof in the case where R(k) has an analytic extensions. We will 
remove this assumption in Section [5] thereby completing the proof. □ 

Equivalence of the formula for <5(fco) given in the previous theorem with the one 
given in the introduction follows after a simple integration by parts. 

Remark 5.5. Formally the equation (]5 . 28[) for q can be obtained by differentiating 
the equation (|5.27p for Q with respect to x. That this step is admissible could be 
shown as in Deift and Zhou [llj . however our approach avoids this step. 



Remark 5.6. Note that Theorem \5.S\ does not require uniform boundedness of the 
associated integral operator in contradistinction to Theorem \A.b\ We only need the 
knowledge of the solution in some small neighborhoods. However it cannot be used 
in the soliton region, because our solution is not of the form I + o(l). 

6. Analytic Approximation 

In this section we want to present the necessary changes in the case where the 
reflection coefficient does not have an analytic extension. The idea is to use an 
analytic approximation and to split the reflection in an analytic part plus a small 
rest. The analytic part will be moved to the complex plane while the rest remains 
on the real axis. This needs to be done in such a way that the rest is of 0(t -1 ) and 
the growth of the analytic part can be controlled by the decay of the phase. 

In the soliton region a straightforward splitting based on the Fourier transform 

(6.1) R(k) = [ e ikx R{x)dx 

JR 

will be sufficient. It is well-known that our decay assumption (|4.32p implies R € 
L 1 (E) and the estimate (cf. [271 Sect. 3.2]) 



(6.2) \R(-2x)\ < const q{r)dr, x > 0, 

J X 

implies x l R(—x) G L 1 (0,oo). 

Lemma 6.1. Suppose R € L 1 (K), x l R(—x) € L 1 (0,oo) and let e,(3 > be given. 
Then we can split the reflection coefficient according to R(k) — R a ^(k) + R Tt t(k) 
such that R a ,t(k) is analytic in < Im(fc) < e and 

(6.3) \R ait {k)e-^\ = 0(t~ l ), < Im(fc) < s, \R rit (k)\ = 0(t~ l ), k e K. 

Proof. We choose R a ^{k) = J^° K r t ^ e lkx R(x)dx with K(t) — ^-t for some positive 
< p. Then, for < lm(k) < s, 



DC 



\R a t{k)e- 0t \ < e' 131 / \R(x)\e- Ira{k)x dx < e- fn e K ® e \\R\\ 1 = ||^|| x e~ C/3 "^° >* , 

which proves the first claim. Similarly, for Im(fc) = 0, 

f°° x l \R{-x)\ ||:c'fl(- a ;)|| £ i ( o |0o) const 
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□ 

To apply this lemma in the soliton region k £ iR + we choose 

(6.4) /3 = min -Re($(fc))>0. 

Im(fc)— £ 

and split R{k) = R a ,t(k) + R r ,t(k) according to Lemma RTT1 to obtain 

(6.5) b±(k) = b at t : ±(k)b r ,t,±(k) = b rtt ,±(k)b a , t ,±(k). 

Here b a ,t,±{k), b rt t,±{k) denote the matrices obtained from b±(k) as defined in (|4.26|) 
by replacing R(k) with R a ^ t (k), R r ^(k), respectively. Now we can move the analytic 
parts into the complex plane as in Section U while leaving the rest on the real axis. 
Hence, rather then (|4.28p , the jump now reads 

iK,t,+(k), k e £+, 

&.,t,-(*) _1 > fce£_, 
b rA -(k)-%t,+(k), k£R. 

By construction we have v(k) = l+0(t~ l ) on the whole contour and the rest follows 
as in Section 0J 

In the similarity region not only b± occur as jump matrices but also B±. These 
matrices B± have at first sight more complicated off diagonal entries, but a closer 
look shows that they have indeed the same form. To remedy this we will rewrite 
them in terms of left rather then right scattering data. For this purpose let us 
use the notation R r (k) = R+(k) for the right and Ri(k) = R-{k) for the left 
reflection coefficient. Moreover, let T r (k, feo) = T(k, ko) be the right and Ti(k, ko) = 
T(k)/T(k,ko) be the left partial transmission coefficient. 

With this notation we have 



(6.7) v(k) = 

where 



L(fc)" 1 6+(fc), Re(fc )<|fc|, 
i?_(fc)- 1 B + (fc), Re(fc )>|fc|, 



K r (-fc)e-""W \ _ II 



b-(k)=r T A -kMr ), b+(k) = 



R r (k), 



T r (k,k ) 2 



and 



B + (k) = -^^M-k)e~^ 
Using (|2.1ip we can further write 

/ i o\ A fl,(fc)c-"t fc ) \ 

(6.8) B_(fc) = ( fl , ( _ fc)c ««w , B + (k)=[ L W . 

\ T,(-fc,fc ) a V \° 1 / 

Now we can proceed as before with B±(k) as with b±(k) by splitting Ri(k) rather 
than R r (k). 

In the similarity region we need to take the small vicinities of the stationary 
phase points into account. Since the phase is cubic near these points, we cannot 
use it to dominate the exponential growth of the analytic part away from the 
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unit circle. Hence we will take the phase as a new variable and use the Fourier 
transform with respect to this new variable. Since this change of coordinates is 
singular near the stationary phase points, there is a price we have to pay, namely, 
requiring additional smoothness for R(k). In this respect note that (|4.32|) implies 
R(k) e C'(R) (cf. [23]). We begin with 

Lemma 6.2. Suppose R(k) G C J (R). Then we can split R(k) according to 

(6.9) R(k) = R a (k) + (k-k Q )(k + k Q )H(k), feeS(fco), 

where Ro(k) is a real rational function in k such that H(k) vanishes at fc , — fco °f 
order three and has a Fourier transform 

with xH{x) integrable. 



Proof. We can construct a rational function, which satisfies /„(— k) = f n {k) for k € 

R, by making the ansatz /„(fc) = ^L+i^ E"=o j\(2k )i ( Q j+ i / j j|;)( fc ~ fc o) 3 (fc+fco) : '- 
Furthermore we can choose ay, [3j € R for j = 1, . . . ,n, such that we can match 
the values of R and its first four derivatives at kg, — fco at f n (k). Thus we will set 
Ro(k) = fi{k) with ao — Re(i?(fco)), (3o = Im(i?(fco)), and so on. Since Ro(k) is 
integrable we infer that H(k) € C 4 (M) and it vanishes together with its first three 
derivatives at fco, — fco- 

Note that <I>(fc)/i = 8(fc 3 — 3k^k) is a polynomial of order three which has a 
maximum at —fco and a minimum at fco- Thus the phase $(fc)/i restricted to 
S(fc ) gives a one to one coordinate transform S(fc ) — > [$(fco)/i> $( — &o)/i] = 
[— 16fcg, 16fcp] and we can hence express -ff(fc) in this new coordinate (setting it 
equal to zero outside this interval). The coordinate transform locally looks like a 
cube root near fco an d — &0j however, due to our assumption that H vanishes there, 
H is still C 2 in this new coordinate and the Fourier transform with respect to this 
new coordinates exists and has the required properties. □ 

Moreover, as in Lemma 16. II we obtain: 

Lemma 6.3. Let H(k) be as in the previous lemma. Then we can split H(k) 
according to H(k) — H a>t (k) + H r , t (k) such that H a ^{k) is analytic in the region 
Re($(fc)) < and 
(6.11) 

|#M(fc)e* (fe) * /2 | - 0(1), Rc($(fc)) < 0,lm(fc) < 0, \H r . t {k)\ = O^ 1 ), fc e M. 

Proof. We choose H a ,t{k) = f™ K{t) H(x)e x *Wdx with K{t) = t/2. Then we can 
conclude as in Lemma EU 



|ffa,*(fc)e* Wt/2 | < H^IMe-^)*^^/ 2 ! < \\H{x)\\i < const 

and 



r-K(t) ^ / ,-K{t) 1 X X 

\H r t(k)\ < / \H(x)\dx < const\ / ~~^dx < const — , < const—. 

J-oo V J-oo * K(tyv 2 t 



□ 
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By construction R a ,t(k) = Ra(k) + (k—ko)(k+ko)H a ,t(k) will satisfy the required 
Lipschitz estimate in a vicinity of the stationary phase points (uniformly in t) and 
all jumps will be I + 0(i _1 ). Hence we can proceed as in Section [SJ 

Appendix A. Singular integral equations 

In this section we show how to transform a meromorphic vector Riemann-Hilbert 
problem with simple poles at ik, —ik, 

m+(fc) = m-(k)v(k), fceE, 
(A.l) 

ReSi K m(fc) = lim m{k) [P2 n ) > Res_i re m(fe) = lim m(k) 



r 



m(—k) = m(k) ^ ^ _ 
lim m(i fc) = (l l) 

k — *oo 

into a singular integral equation. Since we require the symmetry condition for our 
Riemann-Hilbert problem we need to adapt the usual Cauchy kernel to preserve 
this symmetry. Moreover, we keep the single soliton as an inhomogeneous term 
which will play the role of the leading asymptotics in our applications. 

The classical Cauchy-transform of a function / : £ — > C which is square inte- 
grable is the analytic function Cf : C\E — > C given by 

(A.2) Cf(k) = / ^rds, k g C\S. 

Denote the tangential boundary values from both sides (taken possibly in the L 2 - 
sense — see e.g. [H eq. (7.2)]) by C+f respectively C_/. Then it is well-known that 
C + and C- are bounded operators L 2 (T,) — > L 2 (E), which satisfy C+ — C_ = I (see 
e.g. [8]). Moreover, one has the Plemelj-Sokhotsky formula (|29j) 

(A.3) C± = i(i^±I), 

where 

(A.4) H/(fc) = 1/ /^ds, fc G £, 

is the Hilbert transform and / denotes the principal value integral. 

In order to respect the symmetry condition we will restrict our attention to the 
set L 2 s (Yj) of square integrable functions / : S — > C 2 such that 



(A.5) /(-fc) = f{k) 



1 

1 



Clearly this will only be possible if we require our jump data to be symmetric as 
well: 

Hypothesis A.l. Suppose the jump data (£,i>) satisfy the following assumptions: 

(i) S consist of a finite number of smooth oriented finite curves in C which 
intersect at most finitely many times with all intersections being transver- 
sal. 
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(ii) The distance between E and {iy\y > yo} is positive for some yo > and 
±i« £ E. 

(iii) E is invariant under k <—> —k and is oriented such that under the mapping 
k I— > — fc sequences converging from the positive sided to E are mapped to 
sequences converging to the negative side. 

(iv) The jump matrix v is invertible and can be factorized according to v = 
bZ b+ = (I — + w+), where w± = ±(b± — I) satisfy 



(A.6) w ± (-k) = -(° 1 jWfc)(j I 



k e e. 



(v) The jump matrix satisfies 

\\w\\oo = ||u>+||z,oo( E ) + ||W-|U<»(E) < OO, 

( A -7) \\w\\ 2 = ||w+|| L 2 (s) + ||ui_|| L 2 (E) < oo. 

Next we introduce the Cauchy operator 

(A.8) (C/)(fe) = -L f f(s)Q K (s,k) 

acting on vector- valued functions / : E — > C 2 . Here the Cauchy kernel is given by 

/ttilJ- \ /-i 3— \ 

(A.9) n K ( s , k) = •+»*»-* ds = s - fc _j i_ Us, 

\ s— i« s—k/ \ s— fe s—'\k/ 

for some fixed i/s: ^ E. In the case n = oo we set 

(A.IO) fiaoM)=(^ -^) rfS - 

and one easily checks the symmetry property: 

The properties of C are summarized in the next lemma. 

Lemma A. 2. Assume Hvvothesis \A.l\ The Cauchy operator C has the properties, 
that the boundary values C± are bounded operators £ 2 (E) — ► £ 2 (E) which satisfy 

(A.12) C+ - C- = I 

and 

(A.13) (C/)(-W) = (0 *), (C/)(W) = (* 0). 

Furthermore, C restricts to L 2 (E), t/ia£ is 



(A.14) (<?/)(-*) = (Cf)(k) r JJ, fceC\E 

/or / <E i 2 (E) or and by (|A.6|) we a/so have 

(A.15) C±(/« tF )(-A) = C T (/ti;±)(*)r5 J), fceS. 

Proof. Everything follows from (|A.11|) and the fact that C inherits all properties 
from the classical Cauchy operator. □ 
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We have thus obtained a Cauchy transform with the required properties. Fol- 
lowing Section 7 and 8 of [3] , we can solve our Riemann-Hilbert problem using this 
Cauchy operator. 

Introduce the operator C w : L 2 S {Y>) — > L 2 (T,) by 

(A.16) C w f = C + (fw_)+C_(fw + ), /6Lj(S). 

By our hypothesis (|A.7|) C w is also well-defined as operator from Lf(E) — + L 2 (T,) 
and we have 

(A. 17) \\C W \\ L 2^ L 2 < const\\w\\ oo respectively \\C W \\ LT ^ L 2 < const\\w\\ 2 . 

Furthermore recall from Lemma 12.61 that the unique solution corresponding to v = 1 
is given by 

mo(fe) = (/(*;) /(-*;)), 

(A.18) /(fc) = ; , \ , , f 1 + ^±^(2 K )- 1 7 2 e**( iK ) 

v ; JV ' 1 + (2K)- 1 7 2 e t *( IK ) V 

Observe that for 7 = we have f(k) = 1 and for 7 = 00 we have /(fc) = jrj^* . In 
particular, /(fc) is uniformly bounded for all 7 G [0, 00] if |fc — i/t| > s. 
Then we have the next result. 

Theorem A. 3. Assume Hvvothesis \A.ll 

Suppose m solves the Riemann-Hilbert problem (|A.1[) . Then 

(A.19) m(k) = (1 - co)mo(fc) + — / **(«)(«;+(«) +io_(«))fi (t (a, 



E 



ere 



u = m+b , 1 = m_6_ a and cq = f — / u(s)(u;+(s) + u>_(s))f2 re (s, 00) ) 

Here {m)j denotes the j 'th component of a vector. Furthermore, [i solves 
(A.20) (I - C w )(ji(k) - (1 - co)mo(fe)) - (1 - c )a,m (A;) 

Conversely, suppose fx solves 
(A.21) (I - C w )(fi(k) - mo(A)) = C w m (fc), 

and 



'" = ' 2~7ri / ^ s )( w +( s ) + W-(s))O re (s,oo)j 7^ -1, 



then m defined via (|A.19[) . wrat/i (1 — Co) = (1 + So) 1 and /j, = (1 + Co) /i, solves 
the Riemann-Hilbert problem (|A.ip and /1 = m±b± . 

Proof. If m solves (|A.1[) and we set n = m±b~^ , then m satisfies an additive jump 
given by 

m + — m_ = fi(w + + W-). 

Hence, if we denote the left hand side of (|A.19[) by rh, both functions satisfy the 
same additive jump. Furthermore, Hypothesis 13.11 implies that fj, is symmetric and 
hence so is fh. Using (|A.13|) we also see that rh satisfies the same pole conditions 
as mo. In summary, m — rh has no jump and solves (| A. 1() with v = I except for the 
normalization which is given by limfc^oo m(ik) — m(ifc) = (0 0). Hence Lemma l2.6l 
implies m — rh = 0. 
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Moreover, if m is given by (|A.19j) . then (1A.12|) implies 
(A.22) m± = (1 - c )m + C±(p,W-) + C±(pw+) 

= (1 - c )m + C w (jjl) ± pw± 
= (1 - c )m — (I - C„,)ju + /x6 ± . 

From this we conclude that p = m±b~ ± 1 solves (|A.20[) . 
Conversely, if p, solves (|A.21[) . then set 

m(k) = m (k) + ^-7 J p(s)(w + (s) + w_(s))ft c (s, k), 

and the same calculation as in (|A~22|) implies rh± = pb±, which shows that m = 
(1 + Co) _1 m solves the Riemann-Hilbert problem (|A.1[1 . □ 

Remark A. 4. In our case mo(fc) € L^°(S), &wt mo(k) is not square integrable and 
so p 6 £?(E) + Lf{T,) in general. 

Note also that in the special case 7 = we have mo(fc) = (l l) and we can 
choose k as we please, say n = 00 such that cq = Co = in the above theorem. 

Hence we have a formula for the solution of our Riemann-Hilbert problem m(k) 
in terms of mo + (I — C w )~ 1 C w mo and this clearly raises the question of bounded 
invertibility of I — C w as a map from L 2 s (Yj) — > L^(E). This follows from Fredholm 
theory (cf. e.g. [37]): 

Lemma A. 5. Assume Hupothesis \A.l[ The operator I — C w is Fredholm of index 
zero, 

(A.23) ind(I-C w )=0. 

By the Fredholm alternative, it follows that to show the bounded invertibility of 
I — C w we only need to show that ker(I — C w ) = 0. 

We are interested in comparing a Riemann-Hilbert problem for which ||wj||oo 
and 1 1 xt? 1 1 2 is small with the one-soliton problem. For such a situation we have the 
following result: 

Theorem A. 6. Fix a contour E and choose n, 7 = 7*, v depending on some 
parameter t£R such that Hypothesis \A.1\ holds. 
Assume that w l satisfies 

(A.24) Hto'Uoo < p(t) and \\w% < p{t) 

for some function p(t) — > as t — > 00. Then (I — C^tj^ 1 : L 2 (fE) — » L 2 .(fE) exists 
for sufficiently large t and the solution m{k) of the Riemann-Hilbert problem (| A. 1[) 
differs from the one-soliton solution m^fc) only by 0(p(t)) 7 where the error term 
depends on the distance of k to E U {±i«:}. 

Proof. By (|A.17|) we conclude that 

= 0[p{t)) respectively \\C w t\\ Lr ^ L 2 = 0(p(i)) 

Thus, by the Neumann series, we infer that (I — C^*) -1 exists for sufficiently large 
t and 

||(1- C^)- 1 -IWli^li = 0(p(t)). 

Next we observe that 

p* - ml = (I — C^C^ml e Lf(E) 
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implying 

HA*- =0{p(t)) and c* = O(p(i)) 

since ||to ||cx> = 0(1) (note /2 = £i = to ). Consequently c = 0{p{t)) and thus 
m t (k) — m (fc) = 0{p{t)) uniformly in fc as long as it stays a positive distance away 
from £ U {±m}. □ 

Acknowledgments. We want to thank Ira Egorova and Helge Kriiger for helpful 
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